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WEIGHTED DECAY ESTIMATES FOR THE WAVE EQUATION WITH 

RADIALLY SYMMETRIC DATA 

PASCHALIS KARAGEORGIS 


Abstract. We study the homogeneous wave equation with radially symmetric data in 
n > 4 space dimensions. Using some new integral representations for the Riemann operator, 
we estimate the A 00 -norm of the solution. Our results refine those of Kubo 13 El in odd space 
dimensions as well as those of Kubo and Kubota [5] in even space dimensions. However, 
our approach does not really depend on the parity of n. 


1. Introduction 

In this paper, we study the homogeneous wave equation 

dfu 0 — A u 0 = 0 in M n x (0, 00 ) 

Uo(x, 0) =<p(x) (1.1) 

d t u 0 (x, 0) = ^(x) 

with radially symmetric data in n space dimensions. Our goal is to obtain a priori decay 
estimates for the solution when generic assumptions are imposed on the initial data. Such 
estimates have been used by several authors in the existence theory for the semilinear wave 
equation d^u — A u = \u\ p , where p > 1. In a later paper, we are going to utilize our results 
to address the existence of global solutions to the more general nonlinear wave equation with 
potential d'fu — A u = \u\ p -V(x)- u, where p > 1 and V{x) is radially symmetric. 

First, consider the solution uq of EH) when the initial data are such that 

Y \%v(x)\ + Y, AX'MI < £(1 + III)- * 1 '- 1 

H<3 M<2 

for some e > 0 and k > 0. In low space dimensions n = 2, 3, sharp decay estimates for u$ 
were obtained by Kubota |6| and independently by Tsutaya HUGH; see also the earlier work 
of Asakura |T| . In these papers, the assumption of radial symmetry was not needed to control 
the L°°-norm of Uq. To obtain similar decay estimates in higher dimensions n > 4, however, 
one has to consider initial data that are either more regular or else radially symmetric. In 
what follows, we focus on the latter case and allow a generic regularity assumption. Thus, 
our goal is to study the homogeneous problem 

d?u 0 — — — -• d r u 0 

r 

Uo(r, 0) 
d t u 0 (r, 0) 

for a fixed integer n > 4. 
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Before we state our main result, however, let us first introduce some notation. Given an 
integer n > 4, we define the parameters a, m according to the formula 

( \ / l 1 > +r) if n is oc fd n 

(“■“) = ( if „ is even. (1 ' 3) 

Note that m > 1 whenever n > 4 and that the sum a + m = (n — l)/2 is independent of the 
parity of n. Also, we shall frequently use the bracket notation (A) = 1 + |A| for each A G M. 
Now, consider the solution to m when the initial data are such that 
z+i 1 

^A S |^ ) (A)| + ^A S+1 |^ ) (A)| < e\ l ~ m (X) m ~ l ~ k (1.4) 

s=0 s=0 

for some e > 0, k > 0 and some integer 1 < l < m. Note that this condition allows the data 
to be singular at the origin and automatically holds in the case that 

1 +1 1 

£>)* Iv ( ' ) (A)|+£<A)‘ +i 1++)! <£<A)~*. 

s=0 s=0 

Our main result in this paper is the following 

Theorem 1.1. Let n > 4 be an integer and define a, m by (USD- Fix an integer 1 < l < m 
and consider functions ip G C Z+1 (M + ) and fit G C Z (M + ) which are subject to (Q. Then the 
homogeneous equation m admits a unique solution u 0 G that satisfies the following 

estimates when D = (d r ,d t ) and (3 is any multi-index with \/3\ < l. 

(a) When 0 < k < m + a — (n — l)/2, 

| D^u 0 (r, t) | < C 0 £r l ~^\~ m ■ (t — r)~^ (t + r y^~ l+m ~ k . 

(b) When k = m + a = (n — l)/2, 

(t + r)' 


\D 0 uo{r,t)\ < C 0 er 




• (t-r)~' m (t + r) 1 


'I -l-a 


1 + In 


I -l-a 


(t-r). 

(c) When m + a < k < 2(m + a) = n — 1, 

\D' 3 u 0 (r,t)\ < C 0 er l ^^~ m ■ (t — r ) m+a ~ k ~^ (t + r) 

(d) When k = 2 (m + a) = n — 1, 

\D l3 u 0 (r, t) | < C 0 £r l ^~ m ■ (t — r )~ m ~ a ~^\ (j _)_ r yP\- l ~ a . ^ _)_ } n ^ _ r )). 

(e) When k > 2 (m + a) = n — 1, 

\D fi u 0 (r, t) | < C 0 £r l - m ~ m ■ (t - r )- m ~ a - |/J| (t + . 

Besides, the constant C 0 that appears above depends solely on k and n. 


When it comes to the decay rates k < (n — l)/2 of part (a), our method leads to slightly 
sharper conclusions which we establish separately, in Corollary 14.,‘-il Also, when n is odd, our 
conclusions may be further improved due to the strong form of Huygens’ principle. Although 
we shall not bother to prove this explicitly, one has uq G C i+1 (R^_) for odd values of n, while 
the estimate of part (c) holds for any decay rate k > (n — l)/2 whatsoever. 
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Let us briefly describe the role of Theorem 11.11 in the existence theory for the nonlinear 
wave equation d 2 u — A u = \u\ p , where p > 1. Since the nonlinearity \u\ p is a C 1 function, we 
require the same regularity on the initial data. Then the special case l = 1 of our theorem 
asserts that the solution uq of the homogeneous equation is 0(r 1-m ) as r —> 0. To control 
this singularity at the origin, one is forced to impose an upper bound on p, as in Our 

theorem actually refines the estimates of in any space dimension n > 4, however, 

one may avoid such precise estimates by utilizing the method-driven upper bound on p. Our 
main contribution here is that, unlike J5], we are placing no restrictions on the decay rate k 
of the initial data in even space dimensions. 

Similarly, let us consider the nonlinear wave equation d 2 u — A u = \u\ p when p > m. If we 
require the same regularity on the initial data, the special case l = m of our theorem asserts 
that the solution u$ of the homogeneous equation is not singular at the origin. In this case, 
the upper bound on p is redundant, but the estimates of Corollary 14.41 are almost necessary 
for an iteration argument to go through. Although the slightly weaker estimates of S3] do 
suffice, those were only obtained for odd space dimensions. 

The interesting feature in our method is that our approach is the same regardless of the 
parity of n. The plausibility of such an approach is not evident from previous considerations 
which depended on various representations of the Riemann operator for the wave equation. 
Here, some new representations will be established to facilitate this kind of an approach. 


Definition 1.2. Given an integer m > 1, define the mth Legendre polynomial by 

1 d m , 


P m {x) = 

and the mth Tchebyshev polynomial by 

(-1 ) m 


2 m m\ dx m 


T m (x) = 


■ VT 


(X 


dr 


x* 


iy 


(i - x 2 y n - 1/2 . 


(2m — 1)!! ' ~ dx r 

Lemma 1.3 (The Riemann operator). Letting z(\,r,t) be the rational function 

X 2 + r 2 — t 2 


z{X,r,t ) = 


2rA 


we define the Riemann operator L as follows. When n is odd, we set 


[ L f}(r,t) = 


1 




I t~r\ 


X^)Gf{X).p m {z{\r,t))d\ 


2 r (n- 1 )/ 2 

with rn — (n — 3)/2. When n is even, on the other hand, we set 

^2 f t+r f 1 y(n-!)/2j(y) 


\ L f](r,t) = 


27rr (n i )/2 < / |t _ r| J z{X rt) 7(J - z(X,r,t) \/T :r 


da dX 


a- 


V2 


<t-r, 0) rl ^(n-l)/2 


T m (cr) 


27TT’( n LA 


l-i yja - z(X,r,t) VT 


da dX 


a" 


(1.5) 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 


(1.9) 


with m = (n — 2)/2. A solution to the Cauchy problem (11.21) is then provided by the formula 


When -0 G C l 


u 0 (r,t) = [Lfi](r,t) + d t [Lp](r, t). (1.10) 

_) and p G C Z+1 (M + ) for some integer l > 1, this solution belongs to C^(M^_). 
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Proof. Our assertion that dnnp defines a C l function will be established through the proof 
of Theorem o where the derivatives of u$ will be estimated. Thus, we need only establish 
the integral representations m and ED for the Riemann operator. One may derive those 
using the equivalent representations of Lamb [2] or Rammaha ; see also pm. For the sake 
of completeness, however, we shall include their derivation. 

The odd-dimensional representation m is precisely the one that appears in |2|. To prove 
the even-dimensional representation (HD, we manipulate the formula 

1 /', p r a<«-=>/=/(A). r - ( " (A - r - ri) d\ d P 

Jo A Jt 2 ~ P 2 J\p—r\ \A - z{X,r, P ) 2 

which appears as equation (6b) in J|. Here, m = (n — 2)/2 is an integer and T m is the mth 
Tchebyshev polynomial (USD. Switching the order of integration, one arrives at 


[Lf](r,t) = 


nr 


r*t+r 


Lf = 


P 


'jl r f n /2 

1 


\t—r\ J\\-r\ \Jt 2 ~ p 2 

/*max(t-r,0) r\+r 


• A (n “ 2)/2 /(A) ■ 


T m (z(X,r, p)) 

V 1 ~ p ) 2 


dpdX 


+ 


P 


irr 


n/2 


|A—r| \Jt 2 ~ p 2 


■ A (ri_2)/2 /(A) • 


T m (z(\,r, p)) 


V 1 - ^( A ’ r ’ py 


dp dX. 


Note that z( A, r, | A ± r|) = =pl, d p z( A, r, p) = —p/ (rA) and also 

t 2 — p 2 = 2rX ■ ^ z(X , r, p) — z( A, r, t)^j. 


Once we now use the substitution a = z(X,r,p ) in the integrals above, we obtain (11.91) . ■ 

In the remaining of this paper, we shall proceed as follows. In section [21 we establish some 
basic facts about the various functions that appear in the previous lemma. In section [21 we 
combine these facts to estimate the Riemann operator and its derivatives. Finally, section 0 
is devoted to the proof of Theorem 11.11 and its improved version, Corollary 14.31 


2. Basic Estimates and Facts 

Our main goal in this section is to collect a few basic facts about the Riemann operator 
of Lemma PI Its integral representation involves the rational function (ED regardless of 
the parity of n, so we intend to focus on this function first. In our next lemma, we estimate 
its derivatives using a rather painstaking approach. However, we do need such an approach 
in order to gain certain cancellations in our subsequent treatise of the Riemann operator. 

Lemma 2.1. With ( r,t ) G arbitrary and z = z(X,r,t ) as in (11.71) . the sharp estimate 

fel < C(i) ■ A2 + r ((y rZ| , A > 0 (2.1) 

holds for each integer i > 0, and the general estimate 

\d i x did?z\<C(i,j 1 k)-X- 1 - i r- 1 - j -(t + r) 2 - k , 0 < X < t + r (2.2) 

holds for all integers i,j,k> 0. In the case that t > r, one also has the auxiliary estimate 

IflA-l < cm ■ 0 < A < 1 + r (2.3) 

with D = (d r ,dt) and (3 an arbitrary multi-index. 
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Proof. Our first assertion m follows trivially once we explicitly compute 


A 2 -t 2 + r 2 


z = 


; d\z = 


A 2 + t 2 — r 2 


; d\z = c(i) ■ 


t 2 — r 2 


if i > 2. 


2rA ’ 2rA 2 ’ A w r\ i+l 

Let us now turn to our second assertion El- According to Leibniz’ rule, we have 

* 3 

\d\d j r d k z\ A ” 1_il r ” 1_il ' I d^ h d j - jl d k ;{\ 2 -t 2 + r 2 ) |. 

*1=0 J1=0 

Since A, r, t < t + r by assumption, we may thus deduce the desired estimate 

* 3 

\d\d J r d k z\ < C + r) 2 - (i - h) - u - jl) - k 

*1=0 ji=0 

< CA -1- * r ~ l ~ j ■ (t + r) 2 ~ k . 

Finally, we prove our last assertion EJ- In the case that f3 = (j, 0), the inequality 

r 2 + 1 1 2 - A 2 1 


(2.4) 


\D p z\ = \& r z\ < C(j) ■ 


\ r 3+ l 


A > 0 


follows by EU> because z remains unchanged when the roles of A and r are interchanged. 
To settle the remaining case (3 = (j, k ) with k > 1, we resort to our general estimate 112.21) . 
Since t > r by assumption, t + r is equivalent to t, so we get 

\D p z\ = \d J r d k z\ < CX- 1 r~ l ~ j ■ t 2 ~ k < CA" 1 r~ j ~ k ■ t 

because t > r and k > 1. This also completes the proof of our last assertion El. ■ 

Corollary 2.2. With 2A* = ( d\ , d r , d t ) and a any multi-index, one has 

\D:z(\,r,t)\<C(a)-\-W 
whenever 0 < t < 2r and \t — r\<X<t + r. 

Proof. In the case that a = ( i,j , k) with j + k > 1, our general estimate (12.21) gives 

\D?z(\,r,t)\ < CX~ 1_i r~ l ~i ■ (t + r) 2 ~ k < CA” 1 ^ < CX^^ 

because r < t + r < 3r and A < t + r < 3r by assumption. 

In the case that a = (i, 0, 0), on the other hand, our sharp estimate (12.11) gives 


\D*z(X,r,t)\ = \d\z(X,r,t)\ < C(i) ■ 


A 2 + |t 2 


< C(i) • 


A + t + r 


< Cx- 


rX i+1 ~ “ rX i 
since \t — r| < A and A + t + r < 2(t + r) < 6r. In either case then, the result follows. I 

In our next corollary, we concern ourselves with the derivatives of l/(d\z). This function 
will arise as soon as we integrate by parts the integrals of Lemma II.Ml namely the ones that 
appear in the explicit representation of the Riemann operator. To estimate its derivatives 
in a rather precise manner, we first recall El and El, according to which 


(M = 


2rA 2 


A 2 + t 2 — r 2 


1^1 < C(i) ■ 


X 2 + 1 1 2 - r 2 
rX i+1 
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Keeping such expressions intact, one does not have to estimate them separately while dealing 
with their product. This is also the main idea in the following proof, where special emphasis 
is laid on the derivatives with respect to A. 


Corollary 2.3. With z = z(X,r,t ) as in (11.71) . one has 

X 2 ~ l r l-j p-l 
(t — r) 1+ i +k 


\%d>%(dxz)- 1 \<C(i,j,k)- 


(2.5) 


whenever 0 < r < t and 0 < A <t + r. 


Proof. Let us set D = ( d r , d t ) and f3 = (j, k) for ease of notation. By repeated applications 
of the chain rule, one obtains an identity of the form 

\P\ i 

D^diz)- 1 = Y, (SazT 1 -' C(/3 1 .....A,i)-n( B ' J ’ a ^). (2.6) 

1=1 Pi -{-...+/3i=j3 q =1 

where each (3 q is a multi-index. Further differentiating with respect to A, one then obtains 


cEElsr*^)- 1 -'! E 

p!+...+Pl=P 


1=1 s =0 


diU(D^z) 

q =i 


(2.7) 


by Leibniz’ rule. Once we write f3 q = (j qi k q ) for each q , the innermost sum becomes 


E 

Pi+...+Pi=p 


diU(D^) 

<7=1 

Meanwhile, our general estimate (1221 ) ensures that 

< C\- 2 ~ Sri ■ (t + r) 2 ~ k ' 


< E E E c W\dT'^apz 

jl+—+jl=j k 1 +...+ki=k si + ...+Si=s q=1 


d s x q+1 d^dt q z 


0 < A <t + r 


so we may combine the last two equations to arrive at 


E 

pi+...+Pi=p 


< CX~ 2l ~ s r~ l - j ■ (t + r) 


2 l-k 


diU(^d x z) 

<7=1 

Next, we insert this fact in (E2D- Since D = (d r , d t ) and (3 = (j, k) by above, we find 

j+k i 

\d\dld k {d x z)- l \ < CJ2J2 l a r S (M _1 "1 • \- 2l - s r~ l - j • (t + r) 2l ~ k . (2.8) 

7=1 s=0 

To handle the latter derivatives, we will employ our sharp estimate m instead of (El- 
First, we use repeated applications of the chain rule to establish the inequality 


|ar s (M _1 “1 <cY,\dxz\- x 


— 1—l—m 


yag + l 


m= 1 


E lift *i 

ai+...+am=i—s q =1 

in analogy with ESI)- When it comes to the product, EB and El combine to give 

'A 2 + t 2 - r 2 ' 


n \ d T i A son 

<7=1 <7=1 


rA a< < +2 


= 2 m C 1 \d x z\ m 

<7=1 
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so we get 

d x z )“ 1 -'| < C\d x z\- l ~ l ■ \ s ~\ 
Inserting this fact in dUHD , we now arrive at 


j+k 

d\did k t (d x z)-^ < CJ2 \d\z\~ 1 - 1 ■ A _2Z_i r~ l ~ j ■ (t + r) 2l ~ k . 
1=1 


According to our computation Q, we also have 


\d\z\ 1 1 


2A 2 r 


A 2 + £ 2 — r 2 


i+1 

< CA a+2 r l+1 




whence 


j+k 

\d\dld k t {d x z)- l \ <CJ2 A 2 "V 1 " J ' • (t + - r) - * -1 . 

i=i 


Since l < j + k within the last sum, this trivially gives 


^idid^dxz)-^ < CA 2_i r 1 ~ j ■ (t + r) j ~\t - r)- l - j ~ k . 

Besides, t + r is equivalent to t whenever t >r, so the desired estimate m follows. ■ 
The last fact we need to treat the Riemann operator in odd dimensions is also the most 
crucial one and appears in our next lemma. Here, we use the integral representation (USD to 
derive some new representations m- The latter are only valid in the interior of the light 
cone, however, they are less singular at r = 0. 


Lemma 2.4. Let n > 5 be an odd integer and a, m be as in (USD- Suppose that f 0 e C^(M + ) 
and fi G C ,+1 (M + ) for some integer 1 < l < m. When t > r, the Riemann operator (USD is 
then subject to an identity of the form 

di[Lfi\(r,t) -£2.r (2.9) 

for i — 0,1 and all integers i < j < l. Here, Pj m denotes a polynomial (ITTTD which vanishes 
with order j at each o/± 1, while Hij denotes a linear operator which acts on functions of X 
and is defined by either (EU2D or (EH- If we let D = (d r , d t ) and f3 be any multi-index, 
then we also have the estimate 




\3 r 3~\P\ f\h\~3 

(t — j-y+j+w 


i+j 

■^A’" +0+ *| / (*)(A)|, 

s=0 


0 < A < t + r (2.10) 


for i — 0,1 and any integer j > 0, provided that t > r as above. 


Proof. First, we shall use induction on j to prove E3 for the case i = 0 and 0 < j < l. 
Explicitly, the polynomials P jrn of interest are given by the formula 


P 


jm 


X) = 


2 m m\ dx m -i 


(x 2 - 1)^ 


US] 


( 2 . 11 ) 
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Note that P 0m is merely the Legendre polynomial (USD- Besides, Pj m vanishes with order j 
at each of ±1, and we also have Pj m = Pj +l m for each j. As for the linear operators H 0 j we 
are going to need, those are given by the formula 




~ 8 _ 1 

d\ d\z(X,r,t) 


(A m +“/(A)), 


( 2 . 12 ) 


where z(X,r,t ) is the rational function (11.71) . We may introduce them for any integer j < l 
and any function / G C Z (M + ). Let us now proceed to the induction argument. Since t > r 
by assumption, the Riemann operator m takes the form 


[Lfo](r, t) 



(n-i)/ 2 p 0in ( z (X, r, t)) dX. 


Moreover, (n — l)/2 = m + a by our definition (11.31) . so the desired identity 

[£/o](M) = LA J\H 0 ,f 0 \(\,r,t) ■ r~ m ~° P, m (z(\ r,t)) d\ (2.13) 

does hold when j = 0. Suppose it holds for some j > 0. As we have already remarked, the 
polynomials (HTTP are such that 


P jm (z(X,r,t)) = Pj +l m (z(X, r, t)) = ^ ^ ■ d x P j+1 , m (z(X, r, t)). 

Integrating (Em by parts, we thus get to replace the leftmost factor in the integrand by 


r ^ ' [#oj/o](A,r,^ = —[Ho,j+ifo\(X, r, t) 

and the polynomial Pj m by Pj + i ;m . This allows us to finish the inductive proof of (12.131) as 
long as no boundary terms arise in the process. In fact, Pj+i, m vanishes at each of ±1 and 
our definition O) gives z[t ± r,r,t) = ±1, hence no boundary terms arise, indeed. 

Next, we establish E3 for the case i = 1 and 1 < j < l. The linear operators H\j we are 
going to need are given by the formula 

[Hijf](X,r,t) = d t [H 0 j f](X,r,t ) - dx^dxz )' 1 ■ d t z ■ [P 0 j/](A,r,f)j. (2.14) 

We may introduce them for any integer j < l and any function / G C Z " I " 1 (M + ). Let us now 
employ the identity (IXT31) we just proved to write 

[Lfi](r,t) = ^ [ t+r [Hojfi](X, r, t) ■ r~ m ~ a P jm (z(X, r, t)) dX. 

^ Jt-r 



Since j > 1, we have Pj m (z(t ± r, r, t )) = P jm (± 1) = 0 by above, so we find that 

{-iy f t+r 


dt[Lfi\(r,t) = 


+ 


-1) 


a t [P 0 i/i](A,r,t) • r m a P jm (z(X, r , t)) dX 

t—r 

j r t + r 

/ d t z • [JSoj/i] (A, r, t) • A, r, t)) dX. 

't-r 
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Recalling our definition dm we may then integrate the latter integral by parts to get 
dt[Lfi](r,t) = f t+r [Hijfi](X,r, t) ■ r~ m - a P jm {z{\,r,t)) dX. 

^ Jt-r 

This is precisely the desired identity HUD for the case i — 1 and 1 < j < l. 

Finally, we turn our attention to the estimate Hui. which asserts that 


\j r j~jo~ko fjo+ko-j 

< COo, ko,j) ■ (t _ r) , +J+Je+t0 ■ E A m+ “ +i |/ W (A)| (2.15) 

for i = 0,1 and each 0 < A < t + r. Allowing the remaining indices to be all arbitrary, our 
plan is to establish the more general estimate 


j 0 ,k„,]) 


\J~ l 0 r J-J0 po-3 

— r y+j+jo+ko 


i+io+j 

E 

8 =0 


m+a+s 


l/ (s) (A)l 


(2.16) 


for i — 0,1 and each 0 < A < t + r. When i 0 = 0, this actually improves (ITT51) by an extra 
factor r k °t~ k °, which is less than 1 since r < t by assumption. In what follows, we may thus 
focus on the derivation of dzunp , instead. 

To prove dump for the case i — 0, we need to check that 


\j-io r j-jo +jo-j EE 

p?yaf°[g 0j /]| <C(to,Jo,*»,!)■ (f _ r)J+J , +tQ . gA m+ “ +, |/ (,| (A)| (2.17) 

for each 0 < A <t + r. Let us proceed using induction on j. When j = 0 and j 0 + k 0 > 1, 
our task is trivial since H 00 f = X m+a f(X ) depends only on A. When j = 0 and j 0 — k 0 — 0, 
on the other hand, an application of Leibniz’ rule gives 


ip 

\d\°[H 00 f]\ < C(i 0 ) ■ Y, X m+a+s ~ i0 |/ Cs) (A)|- 

s=0 


This proves hud when j = 0, so suppose the same estimate holds for some j > 0. In view 
of our definition (HEP , we may then write 


[H 0J+1 f](X, r ,t) = d x ((d x z) 1 • [H 0j f](X,r,t)y 

The derivatives of (<9 A z) -1 were treated in Corollary 12.hi while those of H 0 jf are subject to 
our induction hypothesis dnzp. Resorting to Leibniz’ rule, one may easily use these facts 
to finish the inductive proof of (Hup , so we shall omit the details. 

To establish (Him for the case i — 1, we need to check that 




po -] 

(t — j \ 1 +3+jo+ko 


*0+1 +j 

Y X m+a+s \f {s) (X)\ 

s =0 


for each 0 < A <t + r. Here, we may proceed directly starting with the definition 
[//ij/ll'A.r.fl = -dx^dxzY 1 ■ d,z ■ [//u+'A,?',<ij 


(2.18) 
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we introduced in (12.1411 . The derivatives of H 0 jf are subject to (12.1711 . those of (<9\£) _1 were 
treated in Corollary 12.41 while our general estimate (12.211 applies for the derivatives of z. In 
view of these facts, our last assertion d2HHD is now easy to deduce using Leibniz’ rule. ■ 
Our next step is to establish an even-dimensional analogue of the previous lemma. Here, 
our task is similar but harder, as the representation (USD for the Riemann operator is more 
subtle. Consider, for instance, the function 

U 0m (X,r,t)=[ [a - £(A,r,t)]” 1/2 • -?== da, \t-r\<X<t + r 

J z(X,r,t) V 1 & 

that appears in (USD as an inner integral. At the endpoints A = r ± t, this function behaves 
rather nicely and actually attains the value 7 t/\/ 2- At the endpoint A = t — r, however, the 
integral above happens to diverge. Namely, our definition <Q gives z(t — r,r,t) = —1 and 
this creates a singularity that is not integrable near a = — 1. It is worth noting that we only 
have to deal with the troublesome endpoint when t > r. To handle this case, we shall need 
to introduce and study the functions 

U jm (X,r,t)= f [a - z(X,r,t)] J ~ 1/2 • -^=2= da, \t-r\<X<t + r 

J z(X,r,t) V 1 — cr- 

for all integers 0 < j < m. These will play a role similar to that of the polynomials Pj m we 
introduced in odd dimensions iram 


Lemma 2.5. Fix an integer m > 1 and let T m be the polynomial (11.61) . With z = z( A, r, t) 
given by (ED whenever (r,t) G set 


Ujm(X, r, t) 



z y - 1 / 2 


T m (a) 

y/1 — cr 2 


da 


(2.19) 


for each \t — r\ < X < t + r and each integer 0 < j < m. 

(a) At the endpoint X — t + r, the function Uj m has a zero of order j. 

(b) Letting D* = (d\, d r , d t ) and a be any multi-index, one has 


\DfU 0m {X,r,t)\ < C(a, m) 


X 


r — t + X 


1/2+M 


• A” h 


whenever 0 < t < 2r and \t — r\ < X <t + r. 

(c) Assume that t>r. Letting D = (d r ,d t ), one then has 


\D l3 U jm (X, r,t)\ < C(m) ■ 


X r 


1 / 2 - 1 , 


t + rj y/r - t + X 

for each \t — r\ < X < t + r and each multi-index f3 with \/3\ < j < m. 


( 2 . 20 ) 


( 2 . 21 ) 


Proof. For the first two parts, we use the substitution v = (a — z)/(l — z) to write 


Uj m (X, r, t) = (1 - z) J 


- 1 7A- 1 / 2 


_ T m (g ) 

\J\ — v a/1 + a 


du, a = v + (1 — u)z. 
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In view of Ok we have z — 1 when A = t + r, so part (a) is clear. To prove the estimate 
of part (b), we take j — 0 in the last equation and differentiate to find that 


D?U 0m {\,r,t) = 




■ Di 


T m (cr) 
a /1 + a 


du, a — u + (1 — u)z. 


( 2 . 22 ) 


When it comes to the rightmost factor in the integrand, we have 


Di 


Tmicr) 


M 


E E C(a u ..., a„ s) ■ 


d s 


T m (cr) 


n< t 




y/l + orj V ' d<J s \\/l + a, 

' v S= 1 Q;i+...+Q! s =Q! x v ' q= 1 

Since A, r and t are as in Corollary 12.21 an estimate of the form \D° '' q z | < C(a q )- A I" 9 holds 
for z, hence also for a = u + (1 — u)z because 0 < u < 1. In particular, a is bounded and 


jja I T m (o') 


< C(a,m) • (l + a)" 1/2 " H • A~ |Q| 


,\/TTo¬ 
by above. Inserting this fact in (irm we then arrive at 

1 


\D:U 0m (\,r,t)\ <CA-I Q I- 


• (1 + cr) ^ du, a = u + (l — u)z. 


'o y/u(l - u) 

It is easy to check that |z| < 1 for each \t — r\ < A < t + r, whence 

l + <j=l + u + (l — u)z = 1 + z + (1 — z)u >1 + 0 
within the region of integration. Meanwhile, the definition <o of 0 ensures that 

_ (r + t + X)-(r-t + X)^r-t + X 

2rA - A 

for each A > |r — 1\. Once we now combine the last three equations, we obtain If 2. 201) . 

Next, we turn to part (c). Here, we wish to differentiate U ]m a total of \f3\ < j times. As 
the integrand in (12.191) vanishes with order j — 1/2 at the lower limit of integration, we get 


D% m (X,r,t)= tD^-zy-Vi. T m ((r) ^ 

Jz V 1 - or 

For the leftmost factor in the integrand, repeated applications of the chain rule give 


(2.23) 


101 




\D f3 (a-z) j - 1/2 \ < C{/3) ■ E E (a - z) j - s ~ 1/2 ■ Y[\D^ 

S= 1 Pi +...+/ 3 s =/3 q= 1 

Moreover, one has \z\ < 1 whenever \t — r| < A < t + r, so one easily finds 
(u - z ) j - s ~ 1/2 < (1 - z) j ~ s • (u - 0 )- 1/2 < 2 j ~ s • (o- - 0 )" 1/2 
whenever z < a < 1 and s < \(5\ < j. Combining the last two equations, we then obtain 


i+ 


M 


cr — 0 


p'-V 2 


s c E E u-w I/ 2 'n+V' 


(2.24) 

s=l Pi-\-...-\~p s =P 9—1 

Since t > r by assumption, our auxiliary estimate (ES> is applicable here. In fact, we have 
r 2 + \t 2 — A 2 | < r(r + t + A) < 2 r(r + t), \t — r\ < A < t + r 
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because \t — A| < r for such A, hence (12.811 ensures that 

< C(/5 9 ) - = C.( A ) . 1±A , | t -r|<A< t + r. 


(2.25) 


Once we now combine the last equation with (i Tm . we get 

\P\ / . , \ S / 

\D\a - zy-^KCr-W - z)~ 1/2 (J < C(a - z)~ x l 2 ( 

S=1 ' ' ' 


t + r 
A r 


because A < t + r. Inserting this fact in (12.2511 . we thus arrive at 


\D% m (X,r,t)\<C 


t + r 
A r 


(a 


r 1 / 2 • 


(1 


<7 


21 - 1/ 2 . 


\T m (a)\da. (2.26) 


Note that the polynomial factor is bounded, as \z\ < 1 by above. If we also neglect a factor 
of (1 + a) -1 / 2 from the integrand, the resulting integral becomes 

f (cr — z)^ 1 / 2 ■ (1 — a) -1 / 2 dcr = f u~ 1/2 {l — u)~ 1/2 dv = 7r 

Jz Jo 

by means of the substitution u — (a — z)/(l — z). In particular, (12.2611 leads to 


t + r 
A r 


\D% m {X,r,t)\<C 

Besides, the definition (dD of z is such that 

(r + t + A) • (r - t + A) 

2Ar — 


m 


(i + z) 


- 1/2 


t + r 
2A r 


(r - t + A) 


for each A > / — r|, so we may combine the last two equations to finally deduce (12.211) . ■ 

The functions of the previous lemma are related to the first component of the Riemann 
operator in even dimensions (HD. To treat the second component, we shall need to study 
some similar functions that we introduce below (12.271) . Although their definition resembles 
the one we had before (12m we are now interested in the values 0 < A < t — r. For such 
values, the rational function (ED is no longer bounded and a much more delicate approach 
is needed. In fact, it is only here that the Tchebyshev polynomial ED becomes important. 


Lemma 2.6. Fix an integer m > 1 and let T m be the polynomial (11.61) . With z = z(X,r,t) 
given by ED whenever (■ r,t ) G set 

W im {\,r,t)= [ (a — z) l ~ 1/2 ■ -^=== da (2.27) 

J -1 V1 - cr 

for each 0 < A < t — r and each integer 0 < i < m. 

(a) At the e?idpoint A = t — r, the function Wi m agrees with the function Ui m of (USD 
up to any derivative of order i. 

(b) Assuming that 0 < A < t — r, the estimate 

f r A \ »-o—*+i/2 

\W im (\r,t)\<C(m)- — ■-- - J 2,28) 

\t + r t — r — X J 

holds for all integers 0 < i < io < m. 
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Proof. The only difference between the definitions of U im and Wi m is that the lower limit 
of integration is z(X,r,t) for the former and —1 for the latter. By (11.711 . these two quantities 
agree when A = t — r, so part (a) follows easily. To prove the estimate of part (b), we use 
the definition (USD of the Tchebyshev polynomial T m to first write 


W im (\,r,t) = C(m) / O-z) 1 1/2 • 




d io 

da l ° 


d 


m—i o 


da m ~ io 


— (1 — a 

— 7.n v 


2\m— 1/2 


da. 


In the case that > 1; the expression in square brackets vanishes at ±1. This means that 
we may integrate by parts iq times to obtain 


W im (X,r,t) = C(i 0 ,i,m )/ (a - z) 




'-1 


d 


m—i o 


da m ~ io 


— (1 — a 

— 7.0 ' 


2\m— 1/2 


da. 


Here, the expression in square brackets is integrable on [—1,1] whenever i 0 > 0, so we get 

\Wi m (X,r,t)\ < C(—l — -) i_ *o- 1 / 2 

for all integers 0 < i < iq. Since the definition m of z provides the inequality 

^ r+t+X t—r—X t+r t—r—X 

— 1 — z = --- > --- 

2r X 2r X 

for each 0 < A < t — r, we may then combine the last two equations to deduce flZZBD - I 
Corollary 2.7. Assume that m > 1 and t >r. Letting D = (d r ,d t ), one then has 


\D^W jm (X,r,t)\<C(m)- 


t + r 


1/2—101 ym—j+1/2 


(t - r) m “/+101 y/t-r- X 


(2.29) 


for each 0 < A < t — r and each multi-index (3 with \/3\ < j < m. 

Proof. Let us first differentiate fl2~27l) to get 

D^W jm (X,r,t) = f D\a - z)^ 2 ■ T ^L da. 

J -i yl -a 2 

Using the chain rule repeatedly, one may write the leftmost factor in the integrand as 
D^(a-zy-'' 2 = J2 E C(l3 u ...,l3„s,]).(a-zy— 1 > 2 .fl(D /3 ’z). 

s=l (5i J r... J r(3 s =(3 9—1 

Recalling our definition (12371) . one may then combine the last two equations to write 

m 

D^W jm (X,r,t) = E E C((3 1 ,...,f3 s ,s,j) ■ Wj_ Stm (X,r,t) ■ ik^>- 

5—1 /3 9—1 

This trivially gives 

101 

\D /3 W jm (X, r,t)\ < cE E \Wj- Sj m(X,r,t)\ ■ n \D^z\ (2.30) 

5—1 /3i -\-...-\-{3 s ={3 9—1 

and we shall now resort to our auxiliary estimate (EJ. It is easy to check that 


r 2 + |t 2 — A 2 1 > r 2 + r(t + A) > rt, 


0 < A < t — r 
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because t — A > r for such A. In particular, our auxiliary estimate (El gives 

\D^z\<C(8) r2 + t2 ~ A2 < 9MA. ^ + r + A J + r * x 
q Ar 1^9 l+i 

Employing this fact in (12.301) . we thus arrive at 




\D^W jm {X,r,t)\ < Cr~M - ^ |W 7 -_ i . m (A,r,t)| ■ 


S —1 


t + r + A t + r — A 


A 


(2 


Here, let us note that 0 < \/3\ — s < j — s < j < m since \/3\ < j <m by assumption. 
Case 1: When 0 < A < (t — r)/2, it suffices to show that 

i/ 2— 101 X m ~ -i+i/2 


\D^W jm (X,r,t)\ < C(m) ■ 


t + r 


(t — r) rn— i+101+V 2 


(2 


because t — r — A is equivalent to t — r. For the exact same reason, we have 


\Wj- S:m (X,r,t)\ < C(m) ■ 


A 


m-j-\-s+ 1/2 


y t + r t — r 

by Lemma [2.61 with i = j — s and z 0 = m. Inserting this fact in (12.311) . we then get 

+ \ \ Ttl —J + l/2 / I | \ i | \ 

A \ ' ft + r + X t + r — A 


|Z^H+ m (A,r,f)| 


S=1 


t + r t — r 


t + r 


t — r 


Since t + r ± A < 2 (t + r) whenever A <t — r, this also implies 

\ \ m—j+l/2 / , . \ |/3| 

r A \ / t + r \ 


\D^W jm (X,r,t)\ < Cr~ 


= C 


t + r t — r 

m-j+l/2-\p\ 


t + r 


J ~ r 

yn-j+ 1/2 


(t — r) m ~t+ I0I+ 1 / 2 ' 


In view of our assumption that j < m, we may now deduce the desired estimate ESI- 
Case 2: When (t — r)/2 < A < t — r and r <t < 2r, we need only show that 

rMi/2-|/3| 

|0^ Jm (A.r,t)| < (2 

y/t — r — A 

Since t + r ± A < 2(t + r) < 6r under the present assumptions, we easily get 


101 


|Z^W im (A,r,f)| < O-'fl • ]T |W+_ s>m (A,r,f)| • 


S=1 


by means of (CTD . Moreover, A < t — r < r for this case, so we find that 

101 

\D^W jm (X,r,t)\ < CX-W . ^ |Wj- Sl m(A, r,t)\. 


S =1 


Applying Le mm a f2.hi with z = j — s = z 0 > we also have 

r A 


\Wj- 8jm (X,r,t)\ < C(m) 


t + r t — r — A 


1/2 


< 


C A 1 / 2 


y/t — r — X 
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so we may simply combine the last two equations to deduce the desired estimate (Mi - 
Case 3: When (t — r)/2 < X < t — r and t > 2r, it suffices to show that 

£.1/2-1/?! 

\D l3 W jm (\,r,t)\ < (2.34) 

V t — r — X 

because A is equivalent to fir. For the exact same reason, (I23TD reduces to 


\D^W jm (\r,t)\ < Cr~W ■ £ r, t)| ■ 

and then Lemma f2 . til applies to give 


S=1 


fir — A 


\D p W jm (X,r,t)\ <Cr~W-J2 


S =1 


t — r — A 


ia-j+s+l/2 


fir — A 


for any integer i 0 with j — s < Iq < m. 

Subcase 3a: If it happens that t — 3r/2 < A, the admissible choice i 0 = j — s yields 


\DPW jm (\,r,t)\<Cr-W- 



r 

t — r — A 


1/2 (^)‘ 


C r i/ 2 -l/®l 

y/t — r — A 


because t + r — A < 5r/2 whenever t — 3r/2 < A. 

Subcase 3b: If it happens that A <t — 3r/2, the admissible choice i 0 = j yields 


\DPW jm (\,r,t)\<Cr-M- 



r 

t — r — A 


s+l/2 


t + r — A 


Cyi/2—1/?| 
y/t — r — X 


because t + r — A < 5(t — r — A) whenever A < t — 3r/2. 

This establishes the desired (EH and also concludes the proof of the corollary. ■ 

We are finally in a position to prove an even-dimensional analogue of Lemma 12.41 


Lemma 2.8. Let n > 4 be an even integer and a,m be as in Q. Suppose /o € C^(M + ) 
and f\ G C i+1 (M + ) for some integer 1 < l < m. Also, assume the singularity condition 


i+l-l 


A s |// s) (A)| =°(\- 2m - 2+s ) as A —> 0 


s=0 


for i — 0,1 and some fixed S > 0. When t > r, one then has the analogue 


r*t-\-r 


dt[Lfi\(r, t) = C 3 / [Hijfi](X,r, t) ■ r - m - a U jm (X,r,t) dX 


't—r 


rt—r 


+ Cj / [H tJ f/(X, r, t) ■ r~ m ~ a W jm (X, r, t) dX 


(2.35) 


(2.36) 


of dZHD for i — 0,1 and all integers i < j < l. Here, the functions Uj m and Wj m are given 
by (gunp and (E23). respectively, while the linear operators Hij are the linear operators that 
we introduced in Lemma \%~1\ 
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Proof. For the case i — 0 and 0 < j < l, we shall use induction on j to verify the desired 


[Lf 0 ](r,t)=C, 




[Hojfo](X,r,t) 


' t—r 


■r m a Uj m (X, r,t) dX 


+ C 3 f r [H 0 j f 0 ](X,r,t)-r- m - a W jm (X,r,t)dX. (2.37) 

Jo 

When j = 0, this follows by our various definitions fl 1.311 . (11,911 . (12.1211 . (12.191) and (12.271) . Let 
us now assume the last equation holds for some j > 0. In view of the definition 

U jm {X,r,t)= [ [a - z(X,r,t)y - 1/2 ■ -j== da 

J z(X,r,t) v 1 & 

we introduced in 42 m it is easy to see that Uj m and U 3+ i >m are related by the formula 


U jm (X, r, t ) 


1 

3 + 1/2 


n (\ a ' d\Uj+i m{X, r,t). 
d\z(X, r , t) 


Moreover, the same formula relates Wj m and W 3+ \, m because of our definition (12.2711 . If we 
use this fact to integrate by parts the integrals of (IT371) . we then get to replace the leftmost 
factor in each integrand by 


J+TyS' k 7)' |iWol(v ’ i) ) = JTTJ 2 ' [ *W°KV' t > 

in view of our definition dSH2D - Since the integration by parts also replaces Uj m by Uj+x >m 
and Wj m by W 3+ i im , we may thus complete the inductive proof of (12.3711 by merely showing 
that no boundary terms arise in the process. Now, the boundary term at A = t + r is zero 
because Uj + i jTn vanishes there by part (a) of Lemma 12.51 At A = t — r, we have two boundary 
terms with opposite signs. In fact, the functions Uj + x >m and Wj + i >m agree at that point by 
part (a) of Lemma 12.61 so these two boundary terms plainly cancel one another. As for the 
boundary term at A = 0, that one is given by 


Cj 

J + 1/2 


[Hojfo\(X,r,t) 
d\z(X,r , t) 


W, 


j+l,m 


(A,r, t) 


J A=0 


Using (Esd. dnnp and (E21 together with our singularity condition (IT351) on /o, one can 
easily check that this expression vanishes as well. Thus, no boundary terms arise, indeed. 

Finally, to establish (IPfll for the case % — 1 and 1 < j < l, one proceeds as in the proof 
of Lemma fO Since only minor modihcations are needed here, we shall omit the details. ■ 


3. The Riemann Operator in High Dimensions 

Using the results of the previous section, we shall now study the Riemann operator, which 
is defined by either m or (HU), according to the parity of n. To estimate its derivatives, 
we divide our analysis into two cases that we treat separately in our next two propositions. 
First, we focus on the interior region t > 2 r and provide an estimate that does not involve 
any boundary terms. For odd values of n, a similar estimate appears in [!2J. The main idea 
we use to handle this region coincides with the one in [Hj, where the wave equation with the 
inverse-square potential is studied. However, the approach in |S] is quite different from ours. 
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Proposition 3.1. Let n > 4 be an integer and define a,m by (USD- Suppose /o G C^(M + ) 
and f\ G C Z+1 (M + ) for some integer 1 < l < m. Also, assume the singularity condition (IT35I) 
for i — 0,1 and some fixed 5 > 0. When D = (d r , d t ) and t > 2 r, the Riemann operator is 
then such that 

i+j 


rt+r 


\D^d\[Lff\{r,t)\ < Cfin) ■ 

C 2 (n) ■ r ] 


— \b\—m—a 


X m-J 


—\u\-m 


t-r ( r-t + X ) 
' ‘t—r y2m 


— - V A 

1 —a / -j 


5 + 1—1 I f( S ) 


\frw\dx 


s=0 

i+j 


(t - r)i +m+a J o (t-r - A) 


- TT^-VA 

— r — a ■ i -—' 


s+1 — i | f(s) 


\frw\dx (3.1) 


s=0 


for i — 0,1 and each max(i, \/3\) < j < l. Moreover, one has C 2 (n) = 0 when n is odd. 

Proof. We divide our analysis into two cases, according to the parity of n. 

Case 1: For odd values of n, an application of Lemma 12.41 allows us to write 

(— 1 ) J ’ b+ r 


d l t [Lfi}(r,t ) = 


[Hijfi\( A, r,t) ■ r m a P jm (z( A, r, t)) dX, i<j<l. 


f t—r 


Recall that P jm is a polynomial with a zero of order j at each of ±1, while z(t ± r, r, t) = ±1 
by (HH). In particular, we may differentiate the last equation \ j3\ < j times without creating 
any boundary terms to obtain 

/ t+r 

Y l ^ 1 (Hi j fill ■ r- m ~ a ~ 1*1 • \D^P jm (z(X,r,t))\dX. (3.2) 

_r 01+02+03 =P 

When it comes to the operators LR,-, the general estimate 

(j\j r i—\0i\ f\0i\-j Lff 


■n i: < 


(t — yfi+i+\0i\ 


'Y X m+a+s |/,p)(A)|, 0 < A < t + r (3.3) 


s=0 


is provided by Le mm a 12.41 Linder our assumption that t > 2r, it actually implies 


i+j 


\D 01 [H ij fi]\ < C\-*- j r j ~ w . Y X m+a+s |/f } (A)|, t-r <\<t + r (3.4) 


s=0 


because A, t ± r and t are all equivalent here. For the exact same reason, flT25j) gives 
\D J z(X,r,t)\ <^(7)-^<^, t-r<X<t + r 

for any multi-index 7 . Moreover, one has \z(X,r,t)\ < 1 for such A, so this easily leads to 

\D^P jm (z(X,r,t))\<Cr-^. 

Using the last equation and (E3D, we may then estimate (El as 

rt+r i+j 


\D s a\\L,U}(r.t)\<Cri 


— \b\—m—a 


! t—r 


X m ~ j ' Y< A s+a_ * | // s) (A) | dX. 


s =0 


Since a = 1 when n is odd, this is precisely the desired estimate m with C 2 (n) = 0. 
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Case 2: For even values of n, an application of Lemma 12.81 allows us to write 


nt-\-r 


dt[Lfi](r, t) = Cj / ■ r~ m - a U jm (X,r,t) dX 


't—r 


ft—r 


+ Cj / [Hijfi\( A, r,t) • r m a W jm ( A, r, t) dX, i <j <1. 

Jo 

By Lemma [2.51 the function U ]m vanishes with order j at the endpoint A — t + r. Also, its 
derivatives up to order j agree with those of W ]m at the endpoint A = t — r by Lemma I2.fi! 
Thus, no boundary terms arise upon \f3\ < j differentiations of the last equation, so we get 

/•t+r 


lD e dHLM\ < C 


tt—r 


Y \D 01 [H ij f i } \ . r - m - a -\fo\ • \D^U jm (X,r,t )| d\ 


(h+/32+/33=/3 
t—r 


+ c f Y \D^ l [Hijfi\ \ ' ■ \D^ 3 W jm (X,r ,t) \ dX 

/3i+/3 2 +/3 3 =/3 

= Ai + A2- (3-5) 

Within the first integral, A is equivalent to t ± r as before, so an estimate of the form 

Cyl/2-lftl 


\D^U jm {X,r,t)\ < 


t — r < X < t + r 


\Jr — t + X ’ 

holds by part (c) of Lemma EH Using this fact along with (18.41) . we then find that 


Ai < f7r ?_ ^ _m_a+1//2 


pt+r 




i+j 


‘-r V^t + X jz tl 




Moreover, we have r < t — r < X within the region of integration, hence also 

nt+r \m.— i 


A 1 < Cr j ~ m 


— \o\—m—a 




t- r y/r-t + X Y 


YX s+a+1/2 ~ l \ti s \X)\dX. 


Since a = 1/2 when n is even, the desired estimate (18.11) is thus satisfied by A ± . 

To treat the second integral A 2 , we recall that t is equivalent to fir because t > 2 r by 
assumption. An immediate consequence of (Ed is then 

n\i r i-\Pi\ ,, 

[Ham < (t _ r)l+y ■ £ Xm+a+ ‘ X Wl, 0 < A < t - r 

and an immediate consequence of Corollary 12.71 is 

n r l/2-|/3 3 | . \m-j+ 1/2 1 

\D^W jm (X,r,t)\ < — 7 -r-— -i= - , , 

I Jm\ - (t _ r )m- J+ l /2 - r - A 

Employing these two facts in (Q, we now get 


0 < A < t — r. 


A2 A 


(j r j-\/3\-m-a+l/2 rt-r y 

(t — r y+j+m+ 1/2 


2m 


i+j 


Vt^r^Tx X- 




s =0 
















WEIGHTED DECAY ESTIMATES FOR THE WAVE EQUATION 


19 


Moreover, A <t — r within the region of integration, so this also implies 


A 2 < 


(J r j-\0\~m-a+l/2 

(t — r y+m+l/2 



X 2m 

yjt — r — A 


i+j 

X s +a+ l/2-*|^)(A)| d\ 

s=0 



Since a = 1/2 when n is even, the desired estimate (13.11) is thus satisfied by A 2 as well. ■ 
In our next proposition, we provide a similar analysis for the exterior region t < 2r. The 
most notable difference from our previous estimate is the presence of some boundary terms. 
A more subtle difference is that the regularity of the initial data is not so crucial anymore. 
Namely, our sharpest conclusions occur for the smallest possible value of j in the following 


Proposition 3.2. When t < 2 r, the assumptions of the previous proposition imply 

rt+r \ rm—\3\ *+•? 


\D l3 di[Lf i \(r,t)\ < C[(n) ■ r~ m ~ a 
C' 2 (n) ■ r 


X r 


\t-r\ (r — t + A) 

—m—a fmax(t-r,0) \ 2m, 1 ' •? 


— a -Y, XS+1 ~ l \ri s) W\dX 

s =0 

\ 2m _ 

■ s+l-i | 


(t - r) m +\h\ J Q (t-r - X) 


.^ A ^ 1 -| / W(A)|dA 


s=0 


i+j-l 


s =0 


+ C[(n) ■ r~ m ~ a A m+a_l/3|+s+1_i |// s) (A)| 


A=|t±r| 


(3.6) 


for i — 0,1 and each max(i, |/3|) < j <1. Moreover, one has C' 2 (ri) = 0 when n is odd. 


Proof. We divide our analysis into three cases. 

Case 1: For odd values of n, the Riemann operator <ra is of the form 

i rt+r 

\Lf,](r,t) = - (3.7) 

" J\t—r\ 


To obtain the desired estimate in this case, we shall not have to distinguish between radial 
and time derivatives, so it is convenient to introduce a multi-index 7 of order |y| = \/3\ + i. 
Once we differentiate the last equation, we then get an identity of the form 


D^Lf^t) = l ^ X m+a f t (X ) ■ D^(r- m - a P m (z{X,r,t))) dX 

Z J\t—r\ V 7 

+ m [A m+a /*(A) • r- m - a P m (z(X,r,t )) 

l7i|=bl-i 


A=|t±r| 


= Bi + B 2 . (3.8) 

Since \t — r\<X<t + r within these terms, our assumption t < 2r makes Corollary 12.21 
applicable. Letting = (d\,d r ,d t ), we thus have the estimate 

\Dfz(X, r,t )| < C(a) ■ A _ '“', \t — r| < A < t + r 

for any multi-index a. Besides, |z(A,r, t)\ < 1 for such A, so we easily get 

\DfP m (z(X, r, f))| < C(a, m ) • A~'“', \t — r\ < X <t + r 


(3.9) 
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by repeated applications of the chain rule. This is actually the only fact we need to control 
the right hand side of m- Due to our assumption that t < 2r, it also implies 

\D* (r~ m - a P m (z(\,r,t)))\ <Cr- m ~ a \t-r\<X<t + r (3.10) 

for any multi-index a, as A < 3r here. Applying the last inequality to the integral term Bi, 
let us now recall that |y| = \/3\ + i to hnd 

pt+r pt+r 

|Bi| < Cr~ m ~ a / A m+ a-h\ . |/.(A)| d\ = C r ~ m ~ a / X m ~ 1/31 • A a_i |/-(A)| dX. 

J\t—r\ J\t—r\ 

Since a = 1 when n is odd, the desired estimate m thus follows for B\. Applying dUDD 
to the boundary terms B 2l we similarly get 


x m + a - S1 


\b 2 \< Cr ~ m ~ a 

|7i|.=|^|+i-l si+s 2 =|7i 

®+DI — i 

= Cr~ m ~ a 'y ^ jA m + a -l /3 l +S2+1 -*|/( S2 )(A)| 

S2=0 


A=|t±r| 


A=|t=tr| 


Moreover, \j3\ < j by assumption, so the desired estimate (EH) follows for B 2 as well. 

Case 2: When n is even and t < r, the Riemann operator EH takes the form 

rt+r 

[Lfi\(r,t)=C’ 0 / X m+a fi(X) ■ r~ m ~ a U 0m (X, r,t) dX (3.11) 

J\t-r\ 

with Uo m given by cm This closely resembles the equation EH we used in the previous 
case, although the factor P m {z) has been replaced by Uo m and the value of a has changed 
due to the change in the parity of n. Differentiating directly as before, we now get 


D^[Lf t }(r,t) = B[ + B' 2 , 

where each B' k denotes the corresponding B *. of EH with U 0m instead of P m (z). 
Since t < 2 r by assumption, part (b) of Lemma 12.51 applies to give 

y \ 1/2+1 a | 


\D:U 0m (X,r,t)\<CX~W- 


t + X 


\t — r\ < X <t + r 


(3.12) 


for any multi-index a. This provides the first fact that we are going to need; the second one 
is the inequality A < r — t + A, which trivially holds for the case t < r under consideration. 
Combining these two facts, we see that E3 remains valid when P m (z) is replaced by Uo m . 
In particular, the argument of Case 1 applies verbatim, except for the part where the exact 
value of a was invoked. Said differently, we need only worry about the integral term 

B[ = C' r Xm+a M X ) • D 7 (r~ m - a U 0m ( x , r , t) \ dX. (3.13) 

J\t—r\ V J 

If we combine EH with our inequality A < r — t + A once again, we get 

CWA-H 

\D°U 0m (X,r,t)\ < — \t-r\<X<t + r 

V'r — t + A 
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for any multi-index a. As long as A < t + r < 3r, we then get 

fj r ~ in—a \ 1 / 2 — | ck | 

\D: (r- m - a U 0m (\r,t)) \ < - = = , \t-r\ < \<t + r 

y/r — t + X 

for any multi-index a. Applying this fact to the integral term eh, we now find 

ft+r \ 1/2—ItI 

\B[\< Cr~ m ~ a / . • X m+a |/,:(A)| dX. 

J\t—r\ Vr -t + A 

Since a = 1/2 when n is even and |y| — \/3\ +i by above, we thus obtain 

rt+r \m-\f3\ 

\B[\ < Cr~ m ~ a / --—— • A 1- * I/A A) I dX. 

1 11 “ J\ t -r\ (■ r-t + xy-« IM 

In particular, we obtain the desired estimate (El for this case as well. 

Case 3: Suppose now that n is even and r < t < 2 r. Although our estimate EH behaves 
nicely at the boundary terms A = r ± t we had before, this is not the case for the boundary 
term A = t — r which may now emerge. One may overcome this difficulty using our approach 
in the previous proposition to avoid the boundary terms altogether. However, that approach 
does not lead to the desired estimate, so we shall need to resort to a combination of the two. 
Proceeding to the details, let us fix some £ e C°°(M) such that 

«*) = {J If S | 5/2 } (3-14) 

and introduce the functions 

C 1 (A) = 1-c( F 4t), &(A) = c(A 7 ). (3.15) 

Then the linearity of the Riemann operator allows us to write 

[Lfi}(r,t) = [L(Ci/i)](r,i) + [L(C 2 fi)](r, t) 

as a sum of two terms to be treated separately. 

In order to treat [L(£i/j)], we proceed as in Case 2 . Since £i vanishes when A < 2 (t — r ) 
by definition, the Riemann operator El takes the form 

rt+r 

[L(Cifi)}(r,t) = C' 0 / A m+a £ 1 (A)/j(A) • r~ m - a U 0m (X, r, t) dX. 

J 2(t—r) 

This closely resembles the equation EH we had in Case 2, although the extra factor £i(A) 
is now present in the integrand. When it comes to the first fact EH we used in Case 2, it is 
provided by Lemma 12.51 whenever t < 2r, so it is still valid. As for the inequality we used, its 
analogue A < 2 (r — t + A) now holds since A > 2 (t — r ) within the region of integration. The 
main question is then whether the extra factor £i(A) will bring any changes to our previous 
approach. Although this function is bounded by our definition EH , its derivatives involve 
powers of t — r we did not have before. Nevertheless, the only place where such derivatives 
may occur is at the endpoint A — t + r, and one has t + r > 3(f — r) whenever t < 2 r. Once 
we now recall our definition EH - EH , we see that all derivatives of £x vanish at t + r. In 
particular, one does not have to deal with such derivatives at all. 
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Next, we focus on [L(Q 2 f,)]. For A sufficiently close to zero, £ 2 (A) = 1 by (13. 14(1 - (13. 151) . so 
the functions Q 2 .fi and /) satisfy the same singularity condition as A —> 0. This means that 
we may employ Lemma 12.81 and argue as in the previous proposition to find 

f t+r 


I £>^[£(<2/.)] I < C 


't—r 


Y • r- m ~ a -W ■ \D*U jm \ dX 


Pi+/32+/33= (3 
t—r 


+ c 


= B'( + b:; 


Y Jlij(Q’fi) ■ r- m ~ a - l^l ■ \D^W jm \ dX 

01+02+03=0 

-2 (3.16) 

for i — 0,1 and each max(i, \ f}\) < j < l. Namely, the condition i < j < l allows us to invoke 
the identity of Lemma [2.81 and the condition \(5\ < j to differentiate that identity \/3\ times 
without introducing any boundary terms. 

Let us note that 3 (t — r) < t + r whenever t < 2 r and temporarily assume the estimate 

i+j 


|£' s, p. J (C 2 / i )]l < 


C X’ 


— r y+j+\0i\ 


■ Y X m+a+S2 |// S2) (A)|, 0 < A < 3 (t-r) (3.17) 


S2=0 


whose proof is given below. For the remaining values A > 3 (t — r), the left hand side is zero 
because ( 2 vanishes for such A by ra - flunp - Employing to , one thus arrives at 


B'l < Cr 


r»3 (t—r) 

Y ■ \D^U jm \ ■ A" H/3l l ^ yn+a+s 2 \f{ S2 \x)\ dX. 


't—r 


01+02+03=0 


S2=0 


Since r is equivalent to t + r whenever t < 2r, part (c) of Lemma ESI ensures that 

ex 1 ' 2 - 1&1 


\D^ 3 Uj m (X, r,t)\ < 


t — r < X < 3{t — r) 


y/r — t + X ’ 

because 3 (t — r) < t + r by above. Once we now combine the last two equations, we get 

r 3(t-r) 


Bf < Cr 


E 


'*- r 01+02+03=0 VT i + A 


r 1 02 1 \m-\0i\-\0s\ yN 

■ ^A S2+ a+i/2^|/f 2) (A)|dA. 


S2=0 


Moreover, r > f — r > A/3 within the region of integration, so this also gives 
S'' < Cr -m- 


rt+r \m-\0\ 

. Y X 82 ^ 1 / 2 -* \f{ S 2 \x)\ dX. 


t-r y/r - t + A 


S2=0 


Since a = 1/2 when n is even, the desired estimate (13.611 is thus satisfied by B". 
Next, we focus on B'Q Noting that r < t + r < 3 r by assumption, we get 


\D*W jm (\,r,t)\ < 


C yn-j+ 1/2 1 

(t - r) m "t +W ' y/t-r-X ’ 


0 < A < t — r 


by means of Corollary 12.71 Together with (13.1711 . this allows us to estimate B £ 


as 


B'{< 


Cr~ 


f*t—r 


(t — r) 1 


i+m 


E 




(t — r) lAI+lthl ./+ _ r _ \ 
01+02+03=0 { ’ VI r A S2=0 


\m+ 1/2 *+•? 

. Y X m+a+S2 \f{ S2 \X)\ dX. 
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Moreover, r t — r) 1 < (t — r) * whenever A < t — r < r, so we find 


B'i< 


Cr~ 


pt—r 


X 


2m 


i+j 


(t - r) m+ \P\ J 0 y/t-r - A 


J2 A S2+ o+ 1 /2-*|// S2 )(A)| dX. 


S 2=0 


Since a = 1/2 when n is even, the desired estimate (13.611 is thus satisfied by B'f as well. 

To finish the proof, it remains to establish (13.1711 . According to Lemma 12.41 we do have 


i Dft (%(c 2 /,)ii < (f _ r)i+j+m ■ E Am+ “ + * 

' ' s=0 


d s 

dA 


(Cz(A)/i(A)) 


for each 0 < A < 3(t — r) because 3 (t — r) < t + r by above. For the case r < t < 2r we are 
presently considering, our remaining assertion (EH will then follow once we show that 


i+j 

E A ‘ 

s=0 


^MA)/,(A)) 


i+j 


<C^A S 2 -|/f 2 ) (A)|, 0 < A<3(t-r). 

S2=0 


(3.18) 


Recalling our definition EH-EED, we easily get this when 2 (t — r) < A < 3 it — r) since 

i+j i+j s i+j 


E A< 

s=0 


d s 

dA 


(Cz(A)/i(A)) 


< E A ^ - r )"~ s ■ i/*°wi < c E A " • irwi 


s=0 S 2=0 


■S2=0 


by Leibniz’ rule. When 0 < A < 2(Z — r), on the other hand, our definition (13.1411 - (13.151) is 
such that all derivatives of (2 vanish at A, whence 


i+j 

E A ‘ 

s=0 


d s 

dA 


(Cz(A)/i(A)) 


i+j i+j 

c(s) / 


E A '-|&(A)M/«(A)|<C^A«.|/f(A)|. 


s=0 


s=0 


In any case whatsoever, the desired estimate EH follows and the proof is complete. 


4. Estimates for the Free Solution 


Using the results of the previous section, we now study the solution u 0 of the homogeneous 
equation El- To prove the estimates that Theorem 11.11 asserts for its derivatives, we shall 
also need two elementary facts that we list separately in our next two lemmas. Their proofs 
follow those of similar results in Dm, for instance, so we are going to omit them. As for 
the parameter a > 0 we introduce below, we are merely interested in the special cases a = 1 
and a = 1/2 that we had before (11.31) . 


Lemma 4.1. Let ( r,t ) G 
(X) b dX 


rt+r 


\t~r\ 


(r - t + 


Li be arbitrary. Assuming that a > 0, one has 

Cr a (t + r) b if b > 

Cr a (t + r) _a ( 1 + In if b = 

V (t-r)J 

Cr a (t + r) a (t — r) a+ if b < 


-a 

-a 


for some constant C depending solely on a and b. 
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Lemma 4.2. Let a > 0 be arbitrary. Assuming that b > 0 and t > r > 0, the integral 

rt ~ r X b (A) c d\ 


Bib, c) = 

satisfies an estimate of the form 


'o 


(t -r - xy- a 


C(t — r) a+h ■ ( t — r) c 

if 

b + c > — 1 

C(t — r) a+b ■ ( t — r) -6-1 • (1 + In (t — r)) 

*f 

b + c = — 1 

C(t — r) a+b ■ (t — r) -6-1 

if 

b + c < -1 


where the constant C depends solely on a, b and c. 


Proof of Theorem ilh Let u 0 (r,t ) be the solution of the homogeneous equation (11,21) 
provided by Lemma 11711 Setting fo = if and fi = ip, we may then write 

i 

D 3 u 0 (r,t) = D' 3 [Ly}{r,t) + D s d t [Lp\(r, t) = y D p d l t [Lfi}{r,t). 

i =0 

This is also the expression we wish to estimate when \j3\ < l. Using the notation above, let 
us now write our assumption m on the initial data simply as 

1 i+l 

EE v+1- ‘ l/i'U)! £ eX ‘~ m • (4.1) 

*=0 s =0 

Here, one may readily check that the singularity condition i FTl holds for i = 0,1. 

Step 1: In the interior region t > 2r, we use Proposition I, ‘hi I with j — l to obtain 


\D^u 0 {r,t)\ < Cy-\P\- m ~ a 


rt+r 


X 


m—l 


1 i-\-l 


+ 


_ r (r - t + A) 
r l-\/3\-m r t ~ r \2 m 


dX 


i =0 s=0 
1 i+l 


(f _ r y+ni+a 


pt—r \'Zm 

/ 7— ■ y y A s+i -* i/f } (A)i dx, 

Jo C t-r-xy-“ n 


i=0 s=0 

where C 2 = 0 for odd values of n. Once we now combine this with m, we get 


l-D^ol < C\£r 


l—\3\—m—a 


rt+r /u m-J-fc dA 


(A) ? 


+ 


C 2 er 


l—\(3\—m ft—r \l-\-m / \ \m—l—k 


X l+m ( X y 


dX 


t-r (r-t + xy~ a (■ t-r) 


Z+m+a 


(t — r — A) 1 " 


= Ai + A2, 


where A2 = 0 for odd values of n. 

In what follows, we only deal with the decay rates 0 < k < 2(m + a) that correspond to 
parts (a) through (c), omitting the similar approach that one needs for parts (d) and (e). 
Since t + r and t — r are equivalent whenever t > 2r, an inequality of the form 

\D p u 0 \ < C£r l ~ m ~ m ■ (t + r ) m ~ l ~ k (4.2) 

implies each of the three desired estimates. One can easily obtain this upper bound for Ai 
since A is equivalent to fir within the region of integration. As for the integral term A 2 , 
we may assume that n is even and thus focus on decay rates 0 < k < 2 m + 1. Write 

A 2 = C 2 £r l ~W- m ■ (t - r )- l - m ~ a ■ B(l + m,m — l — k) 
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using the notation of Lemma 14.21 Since 2 m — k > — 1 by above, the lemma ensures that 

B(l + m, m — l — k) <C(t — r ) l+m+a ■ (t — r ) m ~ l ~ k . (4.3) 

Once we now combine the last two equations, we may deduce the desired (14.211 . 

Step 2: In the exterior region t < 2r, Proposition 13.21 with j = max(i, \ fi\) similarly gives 

rt+r \ i—1/3| /\\m-l-k n 

- / A W “ A + C[er- m ~ a 1 1 ± r\ a+l ~W (t ± r) m ~ l ~ k 
J\t—r\ ( r t + X) a 


\D l3 u 0 (r , t )| < C[er~ 


+ 


C' 2 er~ m ~ a /■ max(t - r ’°) X l+m (X) m ~ l ~ k dX 


(t — r) m+ 1^1 

= ®l + ®2 + ®3, 


(t — r — A) 1 ' 


where B 3 = 0 for odd values of n. Let us only consider the decay rates 0 < k < 2 (m + a), as 
before. In view of HDD, the integral term B3 is then such that 

B 3 < Cer~ m ~ a ■ 1 1 - r\ a+l ~ m u ^ m ~ l ~ k 


(t - r) T 

Since the right hand side appears in the boundary terms B 2 , it thus suffices to treat 

rt+r yl-\p\ t\\m—l—k 


Ei + B 2 = C\er 




w 


dX 


+ \t±r\ a+l ~ m (t±r) 


m—l—k 


\t—r\ ( r-t + X y~ a 

Case 1: If t < 2 r and r < 1, then t + r is bounded and we need only show that 

\D^u 0 (r,t)\ < Cer l -W- m . 

In fact, \t ± r| < 3r for this case, so we easily get the desired estimate 


(4.4) 


rt+r 


(r — t + A) a-i dX + r a 


< Cer l ~ 


Bi + B 2 < Cer l -^- m - a 

\_J \t—r\ 

because \/3\ < l and a > 0 by assumption. 

Case 2: Suppose now that t < 2r and r > 1. Since \/3\ < l by assumption, (14.411 leads to 

rt+r /\\m-k— |, 


Ei + B 2 < Cer 


w 


dX 


—a+( t±r ) 


m-\-a—k— I, 


\t—r\ ( r ~ t + A) ] 

rt+r / \ \ m—k 


< Cer~ m ~ a (; t-r )' 


(A) r 


dX 


t - r I (r-t + X) 


T^ + ( t±r ) 


m+a—k 


and we may use Lemma 14. II to handle the integral. 

Subcase 2a: When 0 < k < m + a, we find that 

B x + B 2 < Cer~ m ~ a (t - r)~ m \r a (t + r) m ~ k + (t± r ) m+a ~ k 


Moreover, m + a — k is positive, so we get 

®i + B 2 < Cer~ m ~ a (t - r)~ W (t + r ) m+a ~ k . 

This does imply the desired estimate, as r is equivalent to (t + r) when r > max(t/2,1). 
Subcase 2b: When k = m + a, our previous approach applies verbatim, although an extra 
logarithmic factor is now included in the estimate that Lemma 14.11 provides. 
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Subcase 2c: When m + a < k < 2(m + a) , Lemma 14.11 yields 

+ B 2 < Cer~ m ~ a (t - r)~ W \r a (t + r)~ a (t - r ) m+a ~ k + (t ± r) m+a_fe 
< C£r- m " a (t - r) m+a - fe -'^l. 

Invoking the equivalence of r with (t + r), we may thus deduce the desired estimate. 

Finally, the uniqueness assertion of our theorem can be established using the same energy 
argument as in [3]. Namely, the condition that 0 imposed was of the form 

\d r u 0 (r,t)\ + \d t u 0 (r,t)\ = 0{r~^ 1 ^ 2+s ) as r —► 0 (4.5) 

for some S > 0. In our case, the estimates we just proved are such that 

\d r u 0 (r,t) \ + \d t u 0 (r,t)\ = as r 0. 

However, l > 1 and m is either (n — 3)/2 or (n — 2)/2, so the condition (14.511 does hold. ■ 


Corollary 4.3. Using the assumptions and notation of Theorem M.il one may improve the 
estimates provided for the decay rates 0 < k < m + a = (n — l)/2 of part (a) as follows. 

(i) When 0 <k<m + a — l, we actually have 

\D 0 u o (r,t)\ < C 0 £r l -W- m ■ (t + r) m ~ l ~ k . 

(ii) When k 0 — 1 < k — m — a + l < k 0 for some integer 1 < k 0 < l, we have 

\D 0 u o (r,t)\ < C 0 er l -W- m ■ (t - r ) mi W-l/3|-fco,o) ^ + ^ m -z-fc-min(H/3|-fco,o) _ 


Proof. If either t > 2r or r < 1, then (t + r) is equivalent to (t — r) and each of the desired 
estimates follows from part (a) of Theorem 11.11 Let us now assume that t < 2r and r > 1. 
As in Case 2 of the previous proof, our task reduces to the estimation of 


= C\er~ 


rt+r X l- 


(A) 


m—l—k 


d\ 


\t—r\ 


(r - t + A) 1 " 


+ \t±r\ a+l ~ m ( t±r) 


m—l—k 


(4.6) 


In what follows, we only concern ourselves with part (ii) because part (i) is easier to settle. 
Given any \t — r\ < A <t + r and any 0 < ko < /, let us write 

X l— \fi\ __ ymin(i— \/3\,ko) . ymax(i— \0\—ko,0) 


Since \/3\ < l by assumption, the exponents on the right hand side are non-negative, so 
yH/ 3 ! < (A) min ( i -l^l -fco ’0) +feo • (t + r) max (H/ 9 l- fe o,°) 

< (\} k ° ■ (t — r ) min l ;_ l /3 l _fc 0.0) _ (3 r )max(H/3|—k 0 ,0) 
whenever \t — r\ < A < t + r and t < 2r. Inserting this fact in (14.611 . we thus arrive at 

]g 1 _)_ ]g 2 < (J £r -m-a+max(l-\f3\-k 0 fi) . U _ ^^min(Z—1/3| — fc 0 ,0) x 
rt+r / \ \m—l— fc+fco J\ 

I \ A ) aA , /+ _l „\m-l-k+k 0 +a 

(r - t + A) 1 - +{t±T> 
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For the decay rates of part (ii), we have m — l — k + k 0 > —a, so Lemma ITT) ensures that 

^ + (t ± r) m -*- fe+fc0+a < C (t + r 
Since r is equivalent to (t + r) when r > max(i/2,1), the desired estimate follows easily. ■ 
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